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Using a controlled analytic non-perturbative treatment, that accounts for the quantum nature 
of the phonons, we derive a model that generically describes cooperative breathing-mode at strong 
electron-phonon interaction in one-band one-dimensional systems. The effective model involves a 
next-nearest-neighbor hopping (that dominates over the nearest-neighbor hopping at strong coupling) 
and a nearest-neighbor repulsion that is significantly enhanced due to incompatibility of neighboring 
dilations/compressions. At non-half filling, upon tuning the electron-phonon coupling, the system 
undergoes a period-doubling second-order quantum phase transition from a Luttinger liquid to a 
conducting commensurate charge-density-wave state: a phenomenon absent in both the Holstein 
model and the t-V model. Using fidelity to study the nature of the quantum phase transition, 
we find that the fidelity susceptibility shows a superextensive power law divergence as well as a 
remarkable scaling behavior: both together establish a second-order transition. 
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I. INTRODUCTION 

Pcrovskitc materials are quite ubiquitous and exhibit 
a variety of interesting and intriguing phenomena such 
as superconductivity or charge ordering (or their co- 
existence), colossal magnetoresistance, ferroelectricity, 
spin-dependent-transport, and the interplay among mag- 
netic, structural, and transport properties^"— . Many ox- 
ides, that have the formula ABO3, assume a perovskite 
structure where two adjacent BOq octahedra share an 
oxygen which leads to cooperative octahedral distor- 
tions. Simple systems that manifest such cooperative 
electron-phonon phenomena are the barium bismuthates 
(BaBiOs). Here, only the 6s electrons are involved in 
transport and these electrons produce only a single nor- 
mal mode distortion, namely, the breathing-mode. In 
pure BaBiOi, the BOq octahedra alternately dilate and 
contract with Bi — O bonds of adjacent octahedra differ- 
ing by about 10% which is indicative of strong electron- 
phonon interaction (EPI)i. Thus the relevant physics is 
dominated by a one-band three-dimensional cooperative 
breathing-mode (CBM). 

There is also compelling evidence of strong EPI 
in manganites (from extended X-ray absorption fine 
structure^ and pulsed neutron diffraction^ measure- 
ments) and in cuprates (through angle-resolved photoe- 
mission spectroscopy^). 

In copper oxides, as pointed out in Refs. [3 and|^, the 
dynamics of the Zhang-Rice singlet^ can be described 
by a one-band system with orbitals centered on copper 
sites. Furthermore, the onsite energy is modulated by the 
movement of oxygen closer or further from the neighbor- 
ing copper ion. Thus the breathing-mode is relevant to 
describe the linear modulation of the onsite energy. Con- 
sequently the copper-oxide planes represent a one-band 
two-dimensional CBM system. 

In the context of the two-band Jahn- Teller manganite 
systems as well, when C-type antiferromagnetism man- 



ifests [as in Lai-^Sr^MnOs for 0.65 < a; < 0.9iS], the 
dz2 orbitals participate in the C-chain ordering. A fer- 
romagnetic C-chain can be looked upon as a one-band 
(i.e., dz2 orbital band) and one-dimensional (ID) CBM 
system that is however Jahn- Teller coupled to neighbor- 
ing C-chains whose spin alignment is antiparallel. 

Understanding the CBM phenomena, in systems such 
as the bismuthates, the cuprates, and the manganites is 
still an open question. The main purpose of this pa- 
per is to study the CBM physics in the simpler case of 
a one-band ID system by taking account of the quan- 
tum phonons [see Fig. [IJb)]. In fact, a controlled an- 
alytic treatment of the many-polaron effects produced 
by quantum phonons in a one-band ID Holstein model 
[see Fig. [Ua)}ii (which is a simpler non-cooperative EPI 
system) has been reported not long agoi^^i^. However, 
definite progress has been made in numerically treating 
the Holstein model at half-filling (by employing a variety 
of techniques)^"— and, to a limited extent, away from 
half-filling^a. 

Owing to its cooperative nature, the EPI leads to non- 
local distortion effects which can change the very nature 
of long range order. While a weak interaction is amenable 
to a Migdal-type of perturbative treatment, the strong 
interaction (even for a one-band system) necessitates a 
non-perturbative approacbi^. As a step towards mod- 
eling CBM distortions in real systems (such as the bis- 
muthates, the cuprates, and the manganites), the present 
work builds up on our previous work on the Holstein 
modeli^ to obtain the effective Hamiltonian for a one- 
band ID CBM system^i"— . Upon inclusion of coopera- 
tive effects in the strong EPI, we show that the system 
changes its dominant transport mechanism from one of 
nearest-neighbor (NN) hopping to that of next-nearest- 
neighbor (NNN) hopping while the effective NN electron- 
electron interaction becomes significantly more repulsive 
due to incompatibility of NN breathing mode distortions. 
Away from half-filling in rings with even number of sites 
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FIG. 1. Molecular chains with d^2 orbital hopping sites (filled 
circles) and oxygen sites (empty circles) in (a) Holstein model, 
(b) one-band CBM system. 



(while the Holstein system without cooperative effects 
remains a Luttinger liquid at all interaction strengths), 
our model (at strong interaction), produces a commen- 
surate charge-density-wave (CDW) state which is sur- 
prisingly conducting and whose period is independent of 
density. Furthermore, using scaling of the fidelity suscep- 
tibility (FS), we demonstrate that the CDW transition is 
a second-order quantum phase transition (QPT). 

The paper is organized as follows. We derive an effec- 
tive polaronic Hamiltonian starting from a CBM model 
in Sec. II. Next, we present the relevant formulae for 
the density-density correlation function and the structure 
factor in Sec. III. We then analyze the strong-coupling 
limiting case of our CBM model in Sec. IV. The nature 
of the QPT and the long range order in our CBM model 
are discussed in Sec. V. Finally, we close in Sec. VI with 
our conclusions. 



II. EFFECTIVE POLARONIC HAMILTONIAN 

To bring out the essential physics, we begin with a ID 
model of spinless electrons hopping in a one-band system 
of dz2 orbitals which are coupled to the oxygens in be- 
tween, via CBM as shown in Fig.[ljb)i. The Hamiltonian 
is expressed as H = Ht + Hep + Hi where the hopping 



term Ht, using standard notation, is given by 



Ht = -t^{c]cj+i + H.c. 



(1) 



with Cj (ri) being the destruction (creation) operator of 

an electron in a d-'a orbital (at site j). The EPI term H^p 
is expressed as 



Hep = -gLL>o\/2Mujoy ] ; 



jih 



(2) 



where g is the electron-phonon coupling (EPC), n^ = 
cjci, qi = Ui — Ui-i represents the expansion of the 
oxygens around the d*2 orbital, and the right-hand-side 



(RHS) oxygen displacement Ui = {al + ai) / y/2M loq. Fur- 
thermore, the lattice term Hi representing simple har- 
monic oscillators is of the form 



u,^-±T,A 



jL^ 3 



^E^' = ^oE«> 



2M 



(3) 



The main difference between the Holstein model and 
the above cooperative Hamiltonian is that in the Hol- 
stein model electrons at different sites are coupled to 
different on-site molecular distortions whereas in the 
present system the electrons on adjacent sites are cou- 
pled to the displacement of the same in-between oxygen. 
Thus in our system, to produce an effective polaronic 
Hamiltonian, we need to devise a modification of the 
usual Lang-Firsov transformational so as to take into ac- 
count the cooperative nature of the distortions. To meet 
this end, we used the following canonical transformation 
H = exp(S')iJcxp(— 5) where S now contains the differ- 
ence in densities on adjacent sites 

S = g^{aj- a] ) (uj -n^+i). (4) 

j 

Then, one obtains II = Hq + Hi where 

Hq =u)q 2, 'A^j ~ 2(7^wo 2_. % + 2g^ajo /^ UjUj+i 
j j j 

-ie-39^^(ctc,+i-fH.c.), (5) 

3 

and 

3 3 

with T± = exjp[±g(2aj — aj_i — aj_)-i)]. On account 
of the cooperative nature of the EPI we obtain an ad- 
ditional term 2g^uio ^ UjUj+i involving NN repulsion 
in Hq and the perturbation Hi now involves phonons 
at three sites as opposed to phonons at only two sites 
as in the non-cooperative case. We consider the case 
texp[— 3(7^] << uJo and perform second-order perturba- 
tion theory similar to that in Refs. [l2 and [2^ The 
eigenstates of Hq relevant for perturbation theory are 
\n,m) = \n)ei ® \m)ph where NN occupied electronic 
states are projected out and |0,0) is the ground state 
(GS) with no phonons. The corresponding eigenenergies 
are En^m = E'f^^ +E^. The treatment to perform second- 
order perturbation theor y is an extension of the method 
followed in Refs. [13 and [2^ and yields the following ef- 
fective Hamiltonian for polarons: 



H 



(2) 



EE' 



\phHu\m)ph{m\phHij\0)ph 
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ph 



E: 



ph 



(6) 



Here (as shown by using Schrieffer- Wolff transformation 
in appendix A of Refs. |25| and 1261). it must be mentioned 
that, when iexp[— 3(7^] << wq, -ffo + H^"^^ represents the 
exact Hamiltonian up to second-order in perturbation 



(even for finite anti-adiabatic values t/ujQ < 1); the small 
parameter [t/ [guiQ)] of the perturbation will be derived 
below. In the above equation ([6]), unlike in Ref. [ij, the 
term Hj produces phonons at sites j, j — 1, and j + 1. 
Hence, we get non-zero contributions only when the index 
i=j — 2, j — 1, j, j + l,orj + 2. Then after some tedious 
algebra one obtains 

^ {[nj(l - nj+i) + (1 - nj)nj+i] 
j 

[F3(4, 1, 1) + 2i^2(4, 1) + i^i(4) + 2F,{1) + F2{1, 1)] 

+ [c]„i(l - 2n,)c,+i + H.c.][2Fi(2) + ^2(2, 2)] 
H.c.]Fi(l) 
H.c.]Fi(-l)}, (7) 



2[Cj_2Cj_iCj^]^Cj 



■j+l^j 



+2[c]-iCj-24+i 
where 



F„(ai,...,a„) = 



rni — 1 



E 



("Iff' 



2\Tni 



(an5 ) 



2\m„ 



^^ 7«i!...m„!(mi + ... + m„) ' 



which for large values of g^ becomes Fn ~ 

exp(5^ Er=i "i)/(ff^ Er=i "0 for ElLi"i > 1- In the 
above Eq. ([7]), the last two terms are a direct consequence 
of the cooperative nature of the EPI and are negligible 
for large g^. More importantly, the relative importance 
of the various coefficients is noticeably different from the 
case where no cooperative effect exists (as explained be- 
low). For large g^, the effective polaronic Hamiltonian 
simplifies to be: 



H. 



eff 



2g'^ojQ 



Sg'^cjo 



Ag^ujQ 



j 
f H.c.) 

^[c]_i(l-2n,)c,+i-HH.c.]. (8) 



-^^-^^T(4^.+i 



t\--'a' 



Notice that the coefficient of the NN hopping is signifi- 
cantly smaller than the coefficient of the NNN hopping 
for large g^ and not-too-small t/uiol This is a key fea- 
ture resulting from cooperative effects. The above effec- 
tive Hamiltonian may be contrasted with the following 
Hamiltonian i?e// for the case where there is no cooper- 
ative EPI [i.e., Hep = —V2guJo J2i '^i(a|+ai)] [see Ref.[l3 
andFig.nja)]: 



H. 



eff 



-2g^uJo > Hi 



3 



2g'^uJo ^ 



nj+i(l-nj) 



-te 



^2^^5:(ctc,+i-fH.c.) 



2„-2s^ 



t^e 



Sff^wo 



5][c]_i(l-2n,)c,+i-|-H.c.]. (9) 



We now provide an explanation of the above results. 
In Eq. dni), the coefficient of the ^^n^j^xiX — rij) term 
can be understood as resulting from a hopping process 
where an electron at site j -I- 1 hops to a neighboring 
site j and back, but the lattice has no time to dis- 
tort (relax) locally at site j (j + 1) and thus yields the 
second-order perturbation energy — i^/(energy change) 
(see Fig. 2 of Ref. 1251 ). On the other hand, the coeffi- 
cient of '^Acj_i{l — 2nj)cj+i + H.c] results when, in 
the intermediate state, site j does not distort/relax dur- 
ing hopping and thus yields texp[— 2g^] x ^ |^ where 

<exp[— 2g^] is due to time I ^.-^ 

site j + 1 (see Fig. 2 of Ref. [f 

cooperative case, the NN hopping dominates over the 

NNN hopping in the small polaron limit. 

Using the above logic we see that the higher order 
terms in perturbation theory, for both cooperative and 
non-cooperative cases, are dominated by the process 
where an electron hops back and forth between the same 
two sites. The dominant term to fcth order is approxi- 
mately given for even k by 



taken to distort the 
In the above non- 
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ffWo 



^?ij+i(l-nj): 



while for odd k by 
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9^0 



k-l 



^(cJcj+i+H.c), 



(10) 



(11) 



where 7 is 2 for the non-cooperative case and 3 for the 
cooperative one. Since each term in the perturbation 
theory should be smaller than ujq, we see that the small 
parameter in our perturbation theory is t/{gujo)^. 

Here, a few observations are in order. Firstly, the coop- 
erative effects, unlike in the Holstein model's case, raise 
the potential of the site next to an occupied site and 
thus make it unfavorable for hopping. Consequently, in 
Eq. (JH) as compared to Eq. ([9|) , the exponent is larger for 
the NN hopping and also the denominators of the coeffi- 
cients arc similarly larger for the hopping-gcnerated NN 
interaction and for the NNN hopping. Next, Lau et alJ- 
obtain the same energy expression for a single polaron 
as that given by Eq. ([5|) (when Uj = 0). Additionally, 
in Ref. 1271 . the authors explain the ferromagnetic insu- 
lating behavior in low-doped manganites by using the 
non-cooperative hopping-generated NN interaction [i.e., 
second term on RHS of Eq. ([9])] after modifying the hop- 
ping term for double-exchange effects. From Eq. ([5]), we 
see that cooperative phenomenon must be taken into ac- 
count as it reduces the ferromagnetism generating inter- 
action st reng th by a factor of 1.5. Lastly, the authors 
of Refs. [2^ and [23 study the formation of bipolarons 
using Frohlich polarons with spin degrees of freedom; al- 
though they use a Lang-Firsov transformation followed 
by a Schrieffcr- Wolff transformation (which is similar to 
our type of perturbation theory), they nevertheless do 



not consider the dominant NNN hopping effects which 
are central to our treatment. 

In the next few sections, we wih analyze the effective 
polaronic Hamiltonian given by Eq. ([5]) and show that 
there is a period-doubling QPT from a Luttinger liquid 
to a conducting CDW when the coupling g increases at 
fixed adiabaticity t/ujQ : the transition is a consequence of 
enhanced NNN hopping and pronounced NN repulsion. 
We employ a modified Lanczos algorithni^S [and use an- 
tiperiodic (periodic) boundary conditions for even (odd) 
number of fermions] to study the QPT in the system. 
In all our numerical calculations involving the effective 
polaronic Hamiltonian, we used the series F„(ai, ...,a„) 
given in Eq. ([7]) and not the approximate coefficients in 
Eq. (HI). 



III. DENSITY-DENSITY CORRELATION 
FUNCTION AND STRUCTURE FACTOR 

In this section, to characterize correlations and analyze 
QPT, we present the relevant formulae for the density- 
density correlation function and the structure factor. The 
two-point correlation function for density fluctuations of 
electrons at a distance I apart is given by 



Nr, 



with filling- fraction (FF) (rij) = -jf- where N is the total 
number of sites and Np is the total number of electrons 
in the system. Then the structure factor, which is the 
Fourier transform of W{1), is given by 



5(fc) = ^e"'^'t¥(/), 



(13) 



where wavevector k 
observe that 



2^ with n=l,2, ,N. Now, we 



^w = fE-E)^W' 

Vieven (odd/ 



with 



and 



E^w 



N 



E^w = - 



'odd 



N 



where Np = Y^- n,- (Nn — Y^, n^) is the number 
operator which gives the total number of electrons at 
even (odd) sites. Hence, we obtain the simple expression 



Sin) 



4((iVe - Nq)^) 

N 



(14) 



We will now analyze the situation where only one sub- 
lattice is occupied and obtain some exact results. When 
we consider odd values of L we note that 



Hence, from Eq. (fT2)) we get 

W^(/odd) = 



0. 



4iV2 
p 



(15) 



Next, we observe that the GS becomes an eigenstate of 
the operators N^ and No with the eigenvalues Np (0) and 
(Np) respectively if the even-site (odd-site) sub-lattice 
is occupied. Consequently, we get 



47V2 
[Sin)] = [5(^)],_ = " 



N ' 



(16) 



where [>5'(7r)]j^j^^ is the maximum value that S{tt) can 
attain. 

To analyze the QPTs, we can treat the rescaled value 
of S{tt) as the order parameter S*{tt) defined as follows: 



S*{n) 



Sin) - [^(Tr)], 



[Sin)] 



[Sin)], 



(17) 



where [^'(Tr)]^^^;^^ is the minimum value of Sin); conse- 
quently, S*in) varies from to 1 during the phase tran- 
sition. 

In the next section, we will study the limiting case 
of large EPC values where the NNN hopping is the only 
relevant transport mechanism in the CBM model leading 
to the t2-V model [see Eq. (|T9l) ]. 



IV. ANALYSIS OF THE t2 - V MODEL - A 
LIMITING CASE OF THE CBM MODEL 

The effective Hamiltonian for the CBM model contains 
three terms, namely, NN hopping, NNN hopping, and NN 
repulsion [as can be seen from Eq. ([5])]. There are two 
possible extreme cases of the CBM model corresponding 
to small and large values of the EPC g. For small values 
of g i^ 1), NN hopping dominates over NNN hopping; 
consequently, Eq. ^ reduces to 

Htv = -t^(c]cj+i+H.c.) + y^njnj+i, (18) 



which is the weU studied t-V model^S*^ with t/V « I 
at the small values of g (~ 1) considered. 

On the other hand, for large values of g, NNN hopping 
dominates over NN hopping and Eq. ([5]) can be simplified 
to 

Ht,v =-t2Xl(4-i(l ~ 2nj)cj+i -fH.c.) 






(19) 



which we shall call as the t2-V model; here, since EPC 
is large (i.e., g > 3), t2/V « 1. However, owing to the 
novelty of the model, we shall study it [i.e., Eq. ([T9| ] 
for arbitrary values of t2/V in rings with even number 
of sites. Next, for t2/V << 1, we note that the system 
always has alternate sites (i.e., one sub- lattice) occupied 
for less than half-filling and above half-filling the other 
sub-lattice gets filled. This can be explained, for less 
than half-filling, as follows. At large repulsion, we shall 
compare the energy for the following two situations: 

1. When there are m^ > (m^ > 0) electrons in 
sub- lattice A (B). 

2. When all the itia + ms = Np electrons are in one 
sub-lattice only. 

In case 1, each electron in sub-lattice B has niB — 1 sites 
blocked in B by other electrons in B and at least [at most] 
niA + l [2?Tiyi] sites blocked in B by electrons in sub- lattice 
A; one can similarly argue for the electrons in sub-lattice 
A. Thus in sub-lattice B(A), each electron can hop to at 

unblocked sites and at least 
1 unblocked sites. 



most Y - mA(B} - mB^A) 



f - '^■rnA(B) - mB{A) 

In case 2, each electron has niA + rriB — 1 sites blocked 
by the other electrons in the same sub-lattice. Hence, 
each electron has exactly -y — ?ti^— to^-I-I unblocked sites 
to hop to. At large repulsion, since case 2 gives electrons 
more number of unblocked sites to hop to, we see that 
the total energy is the lowest when all the electrons are 
present in the same sub-lattice. 

As for the other extreme situation V = 0, for even 
number of electrons, the model has both sub-lattices 
equally occupied. 

In the t2-V model, the ground state energy has a slope 
discontinuity, with the energy increasing up to a critical 
value, after which it is constant for FFs \, ^, and ^ 
[as shown in Fig. [H (a)]. We will now show clearly that 
as the interaction strength increases, at a critical value 
ofV/t2, Ising Z2 symmetry (i.e., both sub-lattices being 
equally populated) is broken and only a single sub-lattice 
is occupied. As depicted in Fig.[2](b), the structure factor 
S{t:) jumps from zero to its maximum value [given by 



Eq. ^ for FFs 



and ^] indicating explosive first- 



order QPT from a Luttinger liquid to a CDW. 

Next, we observe that the number of electrons in even 
and odd sub-lattices are conserved quantities for the t2- 
V model. Therefore, GS of the system is an eigenstate 
of both Ne and No with eigenvalues Ng and No , respec- 
tively. Hence, Eq. p^ simplifies to 



Sin) 



4(JVe - No)' 

N 



(20) 



Then, when Z2 symmetry is respected, for even number 
of electrons Np = 2Ne = 2No, we have S{tt) = and 
for odd value of Np we have ^(Tr) = A/N. We find that 
at a critical interaction strength, as shown in Fig. [31 the 
following dramatic changes occur: (i) the structure factor 
5(7r) jumps from to its maximum value AN^/N; (ii) 






/ 





-2 


/ 


- 


-4 


/ 
/ 


- 


-6 


-"/ .•••■' 


- 


-8 


/,'••■'' 




10 /' 

I.I.I. 


(a)- 



16 


■ Np = 4,N=l6 :~" 


-'—'■- 


14 


■ Np = 6,N=18 ■■■■ ! 


- 


12 


. Np = 8,N=I6 — ■ 1 


- 


10 




- 


8 








6 




- 


4 








2 





- (b). 






2 3 
V 



2 3 
V 



FIG. 2. (Color online) Plots, of (a) ground state energy (E) 
and (b) structure factor value S{Tr) , as a function of interac- 
tion strength (V) for the t2-V model in rings with N sites, Np 
electrons, and hopping t2 = 1. 



W{lodd) also jumps to its large V value of —j^', and 
(iii) W{leven) (for I ^ 0) too jumps and its final value at 
half-filling is 1. For a fixed ^2 and iV, the critical value 
Vq of V increases monotonically as Np decreases. For 
N = 16, Np = 2, and ^2 = 1, we get V^^ k, 4. From finite 
size scaling for half-filling, using Vq — Vp° ex XjN"^ and 
system size N < 20, we obtain Vq^ w 2.83. 

We sec from the above analysis that, at a critical repul- 
sion, the system undergoes a discontinuous transition to 
a conducting commensurate CI? Hastate away from half- 
filling while at half-filling one obtains a Mott insulator. 
Usually commensurate CDW's are insulating (see Ref. 
|32| ) whereas our model surprisingly predicts a conduct- 
ing commensurate CDW. Furthermore, quite unlike the 
Peierls transition, the period of the CDW is independent 
of densityl 



V. ANALYSIS OF THE CBM MODEL 

To analyze the QPT at various FFs of a system gov- 
erned by the effective polaronic Hamiltonian given in 
Eq. ([H, we performed our calculations at values of the 
adiabaticity t/uiQ < 1 and .g > 1 such that the small pa- 
rameter t/{guJo) < 1 and te^^^ ^ ujq. Here we report 
only for the conservative case t/ujQ = 0.1 since the results 
at other values of t/wg < 1 are qualitatively similar (as 
shown in appendix |X| . 

As the value of g increases (in the regime of study 
1 < 5 < 3.5), NN and NNN hoppings compete and the 
system gradually transits from a large-V t-V model to 
a large-V t2-V model; thus, at values of g '--^ 1 wc ex- 
pect the system to be a Luttinger liquid while at g ~ 3 
we should get a CDW. In the next sub-sections we will 
demonstrate that the system indeed undergoes a Lut- 
tinger liquid to a conducting CDW transition with the 
QPT being second-order in nature. The QPT discussed 
in this work is quite different from the metallic Luttinger 
liquid to insulating CDW transition studied by many 
authors^"— in a system with only NN hopping and long- 



range Coulomb interaction. 



A. Study of density-density correlation function, 
structure factor, and order parameter 

First, we calculated W{1) and S{k) at FFs j and i nu- 
merically and the results are displayed in Fig. |4l Upon 
tuning the EPC g, the density-density correlation func- 
tion W{1) gradually changes its nature from decaying to 
oscillatory thereby exhibiting long range order; it then 
attains the value given by Eq. ([T5|) at all odd values of 
I corresponding to the state of only one sub-lattice be- 
ing occupied [see Figs. Slja) and (c)]. Furthermore, the 
structure factor value 5'(7r) increases upon increasing g 
and attains the maximum value given by Eq. ()16p [see 
Figs. |4](b) and (d)]. These observations assert that the 
system undergoes QPT from a Luttinger liquid to a con- 
ducting commensurate CDW state away from half-filling 
with period-doubling. Thus, at a critical value of g, the 
Ising Z2 symmetry (i.e., both sub-lattices being equally 
populated) is broken. Quite surprisingly, our model pre- 
dicts a conducting commensurate CDW without an exci- 
tation gap. Furthermore, similar to the t2-V model, here 
too the period of CDW is independent of density and is, 
in fact, twice the lattice constant. 

We will now compare S(j:) versus g behavior mani- 
fested by our model and the Holstein model in Fig. \E\ 
We see that, while our CBM model appears to undergo a 
QPT, the Holstein model does not seem to do so. We ob- 




FIG. 4. (Color online) Density-density correlation function 
Wil) in the CBM model at ^ = 0.1 and iV = 24 for (a) i- 
filling; and (c) i-fiUing. Structure factor S{k) at (b) i-fiUing; 
and (d) i-filling corresponding to plots of W{1) in (a) and (c) 
respectively. 
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FIG. 3. (Color online) Plots of the density-density correlation 
function W(l) and structure factor S(k) for the t2-V model 
on a N-site ring with Np electrons, interaction strength V, 
and hopping t2 = 1. Structure factor S(k) in (b) and (d) 
correspond to plots of W(l) in (a) and (c), respectively. 



serve that the coefficient of NNN hopping for the Holstein 
model (see Eq. ^) becomes much smaller than that of 
the NN hopping as EPC g increases. Hence, the Holstein 
model, for sufficiently larger values of g, behaves like the 
t-V model; whereas our model can be approximated by 
the t2-V model at large g. We know that the t-V model 
does not undergo a QPT away from half-filling^. There- 
fore, the Holstein model too will not undergo a QPT at 
a non-half FF (which is consistent with the results of 
Ref . I13) . Thus, the Z2 symmetry breaking QPT (at non- 
half filling) in our model is a unique feature which has 
no analog in either the Holstein model or the t-V model. 
Furthermore, at half-filling, the t-V model undergoes a 
QPT when V = 2A^2i^ while the Holstein model suffers 
a QPT at g> l22^; on the other hand our CBM model, 
for the range of EPC g considered (i.e., g > 1), is always 
deep inside the CDW phase since the coefficient of NN 
repulsion is much larger than the hopping terms. 



Plots of the order parameter S* (tt) displayed in Fig. [S] 
also reveal signatures of QPT at FFs j and ^ and at 
different system sizes. Moreover, we observe that the 
increase in S* (tt) becomes sharper as the system size in- 
creases. From the figures it appears that there is either a 
continuous or weakly first-order QPT for both the FFs. 



B. Ground State Fidelity; Fidelity Susceptibility 
and its Scaling Behavior 

Although the order parameter S*{tt) depicts a QPT, 
but the nature of the transition (whether it is first-order, 
second-order, or KT-hke) is not clear. Therefore, we take 
recourse to the study of the ground state fidelity (GSF) 
and FS to characterize the nature of the QPT. The GSF 
is defined as the overlap between GSs at two different but 
near values of the control parameter (say g and g + S) as 
follows: 



F{g,6) = \{M9)\M9 + S))l 



(21) 



where |\['o) 



is the GS of the system and 5 is a small 
2T|) that F{g,S) depends 
on 6. On the other hand, the FS^i, defined below as the 
second derivative of GSF—, 



quantity^. It is clear from Eq, 



XFig) = dlF(g,S)\s=o 



21im 

(5^0 



l-F{g,S) 
<52 



(22) 



is independent of i5. 

The GS of the system, after transition, becomes two- 
fold degenerate as either of the two sub-lattices, namely 
even and odd, can have the larger occupancy. Now, any 
linear superposition of the two degenerate states is also 
a GS. Therefore, the calculated GSF [i.e., the absolute 
value of the overlap of GS I^I^o) at two close by values 
of the control parameter (g and g + S)] becomes arbi- 
trary. To eliminate arbitrariness in the estimate of GSF, 
we start with ^0(5) as our initial guess in the modified 
Lanczos algorithm to get the GS \E'o(5 + 5). 
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FIG. 5. (Color online) Structure factor value S{tv) at — = 0.1 
for (a) i-filling and iV = 16; and (b) i-fiUing and iV = 18 in 
our CBM model and the Holstein model. 




FIG. 6. (Color online) Order parameter S*{-k) in the CBM 
model at ^ = 0.1 for (a) i-fiUing; and (6) i-fiUing. 



Next, we point out a mapping that will enable us to 
perform fidelity calculations in systems with sizes larger 
than the usual sizes accessible to the modified Lanc- 
zos technique. At -^-filling in our CBM model, when 
NN repulsion is much larger than both the NN and the 
NNN hoppings, our model can be reduced to the follow- 
ing model at jr^ — filling but without NN repulsion (for 



N-Np 

similar analyses, see the treatment of the t-V model in 
Ref . ISSi and the mapping of the t-Vi-V2 model in Ref.l26[): 



iJ«7 = -te-3^^5](ctc,+i+H.c.) 



4g'^uJo 



^[c]_i(l-n,)c,+i+H.c.]. (23) 



In the NNN hopping term, because of large NN repul- 
sion, we have ignored the contribution of the sequential 
hopping depicted in Fig. 2(c) of Ref. [23. The above 
prescription reduces the dimension of the Hilbert space 
significantly from ^Cn^ to ^-'^"Cn^- From Eq. ^, we 
also observe that the new effective Hamiltonian contains 
only kinetic terms. Hence, the GS in the CDW phase has 
to be conducting away from half-filling. 

In Fig. [3 we depict F{g,6) and XF{g) as a function 
of g at FFs j and ^. The dip in F{g,5) at the critical 
point increases with the increase in S. This happens be- 
cause of the fact that the distance between two ground 
states in parameter space increases with the increase in S. 
However, XF{g) for different small values of S coincide as 
XF{g) is independent of 5 [as can be seen from Eq. (f22|) ]. 

Additionally, Fig. |5] shows F{g,5 = 0.05), XF{g) and 
Xf„„x (o^ the peak FS) for different system sizes at FFs j 
and i The dip (peak) in F{g, 5) [xF{g)] at the extremum 
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FIG. 7. (Color online) GSF F{g, 5) in the CBM model at ^ = 
0.1 for (a) i-filling and iV = 32; and (c) i-filling and N = 30. 
FS XF{g) for (b) i-filling; and (d) i-fiUing correspond to the 
GSF-plots in (a) and (c) respectively. For the sake of clarity, 
only selected points are shown for 5 — 0.02. 



point increases with the system size TV. Furthermore, for 
a finite system, xf^^^ scales like^^i^ 



XF„ 



N^'. 



(24) 



The logarithmic scale plot of the peak FS value 
XF„,ax(-^) with N shows a linear behavior (see Fig. [He)) 
which confirms a power law divergence of XF„,ax(-^) 
at the extremum point ^max- At large A'^, we obtain 
XF„,^„iN) ~ 7V^™^ at ^-filling; whereas at i-filling we 
gst X-Fmax(A^) ^ 7V^-^^®. The superextensive power law 
divergence of XF,„ax along with the dynamical critical ex- 
ponent value z ^ 1 rule out a KT-likc transition (see 
appendix IB] for details). 

In order to examine the possibility of a second-order 
QPT, we consider the following scaling relation^^i^ for 
Xf{9)- 



{xF^.AN)-XF{g.N)) 
XF{g,N) 



/[A^"(5-gmax)], (25) 



where v is the critical exponent of the correlation length. 
Interestingly, a plot of [xf„„(^) " XF{g,N)]/xF{9,N) 
versus N~[g — gmax), as depicted in Fig. El shows a 
nice scaling relation of xf(3, N) with v taking the values 
1.33 ± 0.01 and 1.41 ± 0.01 for the best fits to the uni- 
versal curves at FFs | and i respectively. The superex- 
tensive power law divergence and the scaling behavior of 
Xf demonstrate that the QPT is second-order in nature. 

Furthermore, as pointed out in Rcfs.l39landl40l average 



FS XF{g)/N around the critical point gc scales like 
XF{g) 1 



TV 



Wc-gr 



(26) 



in the thermodynamic limit, with a being a critical expo- 
nent. The three exponents a, /i, and i/ are related as^i^ 



a 



i^ilJ-l). 



(27) 



The values of the critical exponent a, on using Eq. (|27p . 
turn out to be a ~ 1.33 and a ~ 1.22 for FFs j and 
i respectively. On using finite size scaling, we find the 
critical point gc values to be 2.785 and 2.594 for FFs j 
and g respectively [based on positions of dips (peaks) of 
GSF (FS) in Fig. 0. 



VI. CONCLUSIONS 

We derived an effective Hamiltonian for molecular 
chains involving CBM at strong EPI. The spinless 
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8. (Color online) GSF F(g, S) in the CBM model at 
^ = 0.1 and (5 = 0.05 for (a) i-filling; and (c) i-filling. 
FS XF{g) for (fa) i-filling; and (d) i-filling correspond to the 
GSF-plots in (a) and (c). (e) Plot of the peak values of FS 
XF„ax(N) versus A'^, on a logarithmic scale, at i-filling and 
i-filling and the corresponding power-law fits. 
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FIG. 9. (Color online) Scaling behavior of FS XF{g,N) in 
the CBM model at ^ = 0.1 for (a) ^-filling yielding u = 
1.33 ± 0.01 and for (6) i-filling producing u = 1.41 ± 0.01. 



fcrmion model considered here should be relevant to per- 
ovskite systems with large onsite coulomb repulsion. Our 
analysis shows that our system has an effective Hamilto- 
nian of the form 

Ht-t,-v =-t^(c]cj+i +H.C.) 

3 

-t2^(4-i(l-2n,)c,+i+H.c.) 
i 

+Vj2^j^j+u (28) 



with t2 << t for small g {^ 1), whereas t2 >> t for 
large g (^3); furthermore V is significantly larger than 
both t and ^2 for all values of EPC (1 < g < 3.5) con- 
sidered. Thus, NN and NNN hoppings compete and 
the system transits from a large-V t-V model (with a 
Luttinger liquid GS) to a large-V t2-V model (with a 
period-doubling CDW GS) as g increases; our fidelity 
analysis shows that the QPT is second-order in nature. 
In the past, a density independent charge ordering has 
indeed been observed in manganite systems (see Fig. 2 
in Ref. l4ll). However, since the dimensionality and num- 
ber of bands are different, our findings are not directly 
related to these reported results. Although the reported 
calculations were performed for a conservative value of 
the adiabaticity t/ujo = 0.1, we find that our results are 
qualitatively similar in the whole anti-adiabatic regime 
of t/ujQ < 1 (as shown in appendix A). Furthermore, we 
provide one more model system where the utility of GSF 
and FS in studying the nature of QPT is clearly demon- 
strated. 
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FIG. 10. (Color online) Plots of structure factor value S{n) 
at ^-filling for iV = 16 and -^ = 0.5. 



Appendix A 

In this appendix, we study our CBM model at different 
values of the adiabaticity parameter — and for system 
size iV = 16 and FF |. Firstly, we wish to point out 
that, in the entire anti-adiabatic regime, the Z2 symme- 
try breaking captured by S'(7r) is a novel feature of our 
GBM model which is not present in the Holstein model; 
this claim is endorsed by Figs. [51 and [TUl 

As depicted in Fig. [TTl for various values of ^, the 
order parameter S*{g) rises from to 1 when the con- 
trol parameter g is increased; larger values of — lead to 
QPT occurring at smaller values of g. The GSF F{g, 5) 
and the corresponding FS, at different values of — , are 
portrayed in Figs. H^ a) and [T2T b) respectively. The dip 
(peak) in the GSF (FS) occurs at smaller values of g 
when the adiabaticity — assumes larger values. 

From the above analysis, we find that our model ex- 
hibits qualitatively similar behavior for all values of the 
adiabaticity parameter — < 1. Furthermore, we expect 




FIG. 11. (Color online) Order parameter 5'*(7r) in the CBM 
model at ^-filling, N = 16, and for different values of ^. 
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FIG. 12. (Color online) Plots of (a) GSF F{g,S) at S ^ 
0.05 and (6) corresponding FS in the CBM model at different 
values of ^ and for A'^ = 16 and ^-filling. 



following behavior 



XF„ 



N' 



(B4) 



However, (in contrast to the KT-transition) for a first- 
order or a second-order transition, the numerator of 
Eq. (|B2|) (at the extremum point) does not tend to zero 
as fast as the denominator and hence divergence occurs 
in the thermodynamic limit. Therefore, the divergence in 
Eq. (|B2[) leads to an even stronger power law divergence 
in the FS at the extremum point (as can be seen from 
Eq. mi): 



Xi^„ 



N", 



(B5) 



similar behavior at other FFs and system sizes as well. 



Appendix B 

In this appendix, we will describe a new general ap- 
proach to identify a KT transition. The FS can also be 
written^ as 



y. \{^^{g)\Hj\M9)W 
^0 [Eni9)^Eoigf 



(Bl) 



where Hj is the QPT driving Hamiltonian. On the other 
hand, according to perturbation theory, the second order 
perturbation to the GS energy takes the form 






[Eoig) - E^{g)] 



(B2) 



At the KT-transition point, both the numerator and de- 
nominator of Eq. (jB2p tend to zero with system size in 
exactly the same manner; hence no divergence results. 
At the extremum point, the mass gap typically vanishes 
with system size as 



£'l(ffmax) ~ Eo{g,nax) 



N^ 



(B3) 



where z is the dynamical critical exponent. Thus, for a 
KT-transition the FS at the extremum point exhibits the 
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FIG. 13. (Color online) A log-log plot of the excitation energy 
A(iV) versus A^ in the CBM model at i-fiUing and |-filling 
and the corresponding power-law fits. 



where ^ > z. Now, in Fig. I13[ we depict variation of 
the mass gap A(iV) with iV on a logarithmic scale and 
observe linear behavior. At large N, we obtain A{N) ^ 
j^i\,uv (i.e., z ~ 1.007) at i-filling while at i-filling we get 
A(7V) -- jfTnm: (i-C, z ^ 0.984); whereas from Fig. [DJe), 
we find XF.„JiV) -- iV^°°^ (i.e., /i = 2.001) at ^-filling 



while at i-fiUing we get xf^^A^) ^ iVi-^es (jg.. 



1.868). Clearly, /i > z for both the FFs which rules out 
the possibility of a KT-transition. Hence, the QPT in 
our model is either first-order or second-order. 

Lastly, we would like to mention that another ap- 
proach, involving examining the global geometric entan- 
glement, to detect the elusive KT quantum phase tran- 
sition has been recently rcported^i^. 



^ A. Taraphder, R. Pandit R., H. R. Krishnamurthy, and T. 

V. Ramakrishnan, Int. J. Mod. Phys. B 10, 863 (1996). 
^ T. Hotta, Rep. Prog. Phys. 69, 2061 (2006). 
^ J. van den Brink, and D. I. Khomskii, J. Phys.: Condens. 

Matter 20, 434217 (2008). 
* A. Lanzara, N. L. Saini, M. Brunelli, F. Natali, A. Bian- 

coni, P. G. Radaelli, and S. W. Cheong, Phys. Rev. Lett. 

81, 878 (1998). 



D. Louca, T. Egami, E. L. Brosha, H. Roder, and A. R. 

Bishop, Phys. Rev. B 56, R8475 (1997). 

A. Damascelli, Z. Hussain, and Z. X. Shen, Rev. Mod. 

Phys. 75, 473 (2003). 

For a numerical investigation of single-polaron properties 

of the system in Fig. [Tfb), see B. Lau, M. Berciu, and G. 

A. Sawatzky, Phys. Rev. B 76, 174305 (2007). 

Glen L. Goodvin and Mona Berciu, Phys. Rev. B 78, 



11 



235120 (2008). 

^ F. C. Zhang, and T. M. Rice, Phys. Rev. B 37, 3759 (1988). 

^° A. Szewczyk, M. Gutowska, and B. Dabrowski, Phys. Rev. 
B 72, 224429 (2005). 

" T. Holstein, Ann. Phys. (N.Y.) 8, 343 (1959). 

^^ S. Datta, A. Das, and S. Yarlagadda, Phys. Rev. B 71, 
235118 (2005). 

^^ See also Section V of S. Datta, and S. Yarlagadda, Phys. 
Rev. B 75, 035124 (2007). 

^^ Martin Hohenadler, Holger Fehske, and Fakher F. Assaad 
Phys. Rev. B 83, 115105 (2011). 

^^ J. E. Hirsch, and E. Fradkin, Phys. Rev. B 27, 4302 (1983). 

^® C. E. Creffield, G. Sangiovanni, and M. Capone, Eur. Phys. 
J. B 44, 175 (2005). 

" H. Zheng, D. Feinberg, and M. Avignon, Phys. Rev. B 39, 
9405 (1989). 

^* C. A. Perroni , V. Cataudella, G. De Filippis, G. ladonisi, 
V. Marigliano Ramaglia, and F. Ventriglia, Phys. Rev. B 
67, 214301 (2003). 

1^ R. J. Bursill, R. H. McKenzie, and C. J. Hamer, Phys. 
Rev. Lett. 80, 5607 (1998). 

^° M. Hohenadler, D. Neuber, W. von der Linden, G. Wellein, 
J. Loos, and H. Fehske, Phys. Rev. B 71, 245111 (2005); 
M. Hohenadler, G. Wellein, A. Alvermann, and H. Fehske, 
Physica B 378-380, 64 (2006). 

^^ A different type of cooperative electron-phonon interac- 
tion effect corresponding to a distance dependent electron- 
phonon interaction has been considered for a Frohlich po- 
laron by A. S. Alexandrov and P. E. Kornilovitch, Phys. 
Rev. Lett. 82, 807 (1999). Interestingly, their analytic ap- 
proach (like our method) is also valid for t/ujo < 1. 

^^ The many-polaron problem with infinite-range electron- 
phonon coupling was studied for the spinless Frohlich 
model by A. S. Alexandrov and P. E. Kornilovitch, J. 
Phys.: Condens. Matter 14, 5337 (2002). The cooperative 
effect considered by these authors is long-ranged whereas 
our treatment of the breathing mode involves short-range 
EPL 

^^ For a review of many-polaron effects for both Frohlich po- 
larons and Holstein polarons involving spin and spinless 



fermions, see A. S. Alexandrov and J. T. Devreese, Ad- 
vances in Polaron Physics (Springer, Berlin 2009). How- 
ever, this book does not cover the dominant NNN hopping 
due to cooperative EPI effects considered in our paper. 

^^ LG. Lang, and Yu.A. Firsov, Zh. Eksp. Teor. Fiz. 43, 1843 
(1962) [Sov. Phys. JETP 16, 1301 (1963)]. 

'^^ S. Reja, S. Yarlagadda, and P. B. Littlewood, Phys. Rev. 
B 84, 085127 (2011). 

^^ S. Reja, S. Yarlagadda, and P. B. Littlewood, Phys. Rev. 
B 86, 045110 (2012). 

^■^ G. V. Pai, S. R. Hassan, H. R. Krishnamurthy, and T. V. 
Ramakrishnan, Europhys. Lett. 64, 696 (2003). 

^* A. S. Alexandrov, Europhys. Lett. 95, 27004 (2011). 

^^ A. Alexandrov and J. Ranninger, Phys. Rev. B 23, 1796 
(1981). 

^" E. R. Gagliano, E. Dagotto, A. Moreo, and F. C. Alcaraz, 
Phys. Rev. B 34, 1677 (1986). 

^^ F. D. M. Haldane, Phys. Rev. Lett. 45, 1358 (1980). 

^^ G. Griiner, Density Waves in Solids, (Addison- Wesley, 
Reading, MA, 1994). 

^■^ Sylvain Capponi, Didier Poilblanc, and Thierry Giamarchi, 
Phys. Rev. B 61, 13410 (2000). 

^* G. Fano, F. Ortolani, A. Parola, and L. Ziosi, Phys. Rev. 
B 60, 15654 (1999). 

^^ Martin Hohenadler, Stefan Wessel, Maria Daghofer, and 
Fakher F. Assaad, Phys. Rev. B 85, 195115 (2012). 

^^ P. Zanardi, and N. Paunkovic, Phys. Rev. E 74, 031123 
(2006). 

^■^ W. L. You, Y. W. Li, and S. J. Gu, Phys. Rev. E 76, 
022101 (2007). 

^^ R. G. Dias, Phys. Rev. B 62, 7791 (2000). 

^3 S. J. Gu, H. M. Kwok, W. Q. Ning, and H. Q. Lin, Phys. 
Rev. B 77, 245109 (2008); ibid. 83, 159905(E) (2011). 

■*° S. J. Gu, Int. J. Mod. Phys. B 24, 4371 (2010). 

"^ G. C. Milward, M. J. Calderon, and P. B. Littlewood, Na- 
ture 433, 607 (2005). 

■*^ Roman Oriis and Tzu-Chieh Wei, Phys. Rev. B 82, 155120 
(2010). 

■*'' Jean-Marie Stephan, Gregoire Misguich, and Fabien Alet, 
Phys. Rev. B 82, 180406 (2010). 



